Electric field driven insulator-to-metal phase transitions 
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We show that strong enough electric fields can trigger nucleation of needle-shaped metallic embryos 
in insulators, even when the metal phase is energetically unfavorable without the field. This general 
phenomenon is due to the gigantic induced dipole moments acquired by the embryos which cause 
sufficient electrostatic energy gain. Nucleation kinetics are exponentially accelerated by the field- 
induced suppression of nucleation barriers. Our theory opens the venue of field driven material 
synthesis. In particular, we briefly discuss synthesis of metallic hydrogen at standard pressure. 



PACS numbers: 05.70.Fh, 64.60.qe, 64.60.Bd, 82.60.Nh 

The topic of metal-insulator transitions has long been 
established [3, 0] and remains active. Application as- 
pects aside, it concentrates mostly on the underlying 
microscopic mechanisms, such as those induced by dis- 
order (the Anderson transition), intra-atomic interac- 
tion (Mott- Hubbard), band crossing, and some others. 
The common feature of these transitions is that they are 
driven by changes in some material parameter: degree of 
disorder, doping concentration, etc. 

Here, we introduce the concept of electric field driven 
insulator-to-metal phase transitions. They start with 
needle-shaped metal embryos forming in an insulator 
when the system is immersed in a strong enough field. 
We argue that, even when the bulk metallic phase as 
such is energetically unfavorable, increasing the field will 
eventually cause the transition to occur. Furthermore, 
assuming a fixed electric field, the final state of the sys- 
tem will be a uniformly metallic needle-shaped body. 

This concept holds regardless of the microscopic mech- 
anism of the transition (densification, crystallization, 
electron solvation, or others [H, Q); it applies equally 
to solids and liquids. For example, it predicts conduc- 
tive needle-shaped crystallites forming in an insulating 
glass under strong enough fields. As another example, 
strong electric fields will trigger nucleation of liquid Si 
(metallic) needle-shaped inclusions in a semiconducting 
Si host, even at temperatures well below melting. A more 
provocative example, briefly discussed in this Letter, is 
the field-induced synthesis of metallic hydrogen under 
standard pressure. 

The fact that symmetry-breaking electric fields can 
dramatically affect nucleation processes was recently re- 
alized Q while studying nucleation of highly conduc- 
tive filaments in chalcogenide glasses of phase change 
memory Q. Other related phenomena can include bias- 
induced metal-insulator transitions in resistive random 
access memory [5( , and dielectric breakdown in thin- film 
devices Q . Another category is non-photochemical laser 
induced nucleation Q- 



The theory herein takes a significant additional step 
by predicting that sufficiently strong fields can trigger 
transitions to states that would not be stable without 
the field; they will remain metastable upon field removal. 
As a possible candidate we mention bias-induced switch- 
ing from insulating to highly conductive states, such as 
in vanadium dioxide (V0 2 ) [8j and chalcogenide glasses 
[J] , where the conductive phase disappears upon field re- 
moval. Other known phenomena could be pertinent, such 
as e. g. dielectric breakdown in thin oxides. However, 
mainstream understanding of the latter refers to a build- 
up of defects produced by stress, eventually forming the 
onset of a percolation path across the oxide jfjj]. These 
types of mechanisms are beyond the present framework. 

We start with a brief introduction to the electric field 
effect in classical nucleation theory (CNT). The free en- 
ergy of a new particle in the presence of the field is, 



F = Act- Qfi + F E . 



(1) 



Here A and f2 are the particle surface area and volume, a 
is the coefficient of surface tension, and n is the chemical 
potential difference between the two phases, taken to be 
positive when the bulk new phase is energetically favor- 
able. Eq. ([I]) does not specify the type of transition. The 
electrostatic term has the form 
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where e is the electric permittivity of the host insulat- 
ing phase and the effect of particle geometry is embodied 
in the depolarizing factor, n. For a sphere, n = 1/3, 
A = AttR 2 , and ft — AttR 3 /3. In zero field, the maximum 
of F(R) from Eq. (fTJ) provides the nucleation barrier 
Wo = 167T0- 3 /(3fi 2 ) and radius i?o = 2er/|/x| with typi- 
cal magnitudes near 1 eV and 1 nm. Maintaining the 
assumption of spherical geometry, the field reduces the 
nucleation barrier and radius according to [Icj . 
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where E = 2{W /{eR i )} 1 / 2 is typically in the range of 
several MV/cm. 

In CNT, /i must be positive in Eq. (|TJ) , implying a 
metastable host phase. Another important assumption 
is that of spherical symmetry. We show next how phase 
transitions are possible even for the case of energetically 
unfavorable bulk new phase (negative fj,) when a strong 
electric field is applied and the constraint of spherical 
shape is relaxed. 

In our concept, the free energy of Eq. (JlJ has two 
degrees of freedom: spherical symmetry is broken and, 
when the field is sufficiently strong, needle-like conduc- 
tive particles aligned with the field become energetically 
favorable. That can be understood by comparing the 
electrostatic energy contribution of a sphere to that of a 
prolate spheroid of height H and radius R, for which the 
depolarizing factor is [9(, 



= (R/H) 2 [\n(2H/R) - 1] = (R/H) 2 L. 



(4) 



Considering particles of equal volume, the electrostatic 
contribution is greater for a prolate spheroid by a huge 
factor of approximately (H/R) 2 ^> 1 (see Fig. [1]). Physi- 
cally, this enhancement is due to a large induced electric 
dipole in the needle-shaped particle. Once created, it 
will act as a lightning rod, concentrating the field and 
triggering further nucleation. 
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FIG. 1: (a) Absolute, normalized value of the electrostatic 
energy contribution, Fe/Wo, to the free energy as a function 
of nucleus aspect ratio; H/R — 1 corresponds to a sphere. 
A value of E/Eo — 0.25 and nucleus volume of 47ri?o/3 were 
used for comparison purposes, (b) Surface plots of the elec- 
tric potential and streamlines of the electric field (shown in 
red) illustrate the greater electrostatic energy reducing effect 
of elongated metallic nuclei versus spheres. An elongated nu- 
cleus concentrates the field at its tips, similar to the lightning 
rod effect, possibly triggering further nucleation events. 

The exact shape of the elongated nucleus is not known, 
but modeling with either spheroidal or cylindrical parti- 
cles leads to differences only in numerical coefficients [3| . 
We opt for the mathematically concise cylinder shape 
with A = 2ttRH, = ttR 2 H, and free energy, 



F cy i — 



W ( 3RH 3R 2 H 



E 2 H ] . 



with the approximation n w (R/H) 2 for H ^> R in Eq. 
(|4]). As illustrated in Fig. [2] the free energy landscape 
exhibits a range of low nucleation barriers at small R. 
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FIG. 2: Free energy landscape of FIMS as a function of nu- 
cleus radius, R/Ro, and height, H/Rq [from Eq. (0 with 
E/Eo = 0.25]. Contour spacings are F/Wo ~ 0.1. Regions 
of positive and negative free energy are separated by the zero 
contour (red). In (a), the new phase is stable in the bulk 
[negative sign in Eq. [5], while in (b) the new phase would be 
unstable [positive sign in Eq. [5] if the electric field were not 
present (note the difference in scale) . The contours show that 
when the field is strong enough, nucleation pathways with 
much lower barriers become available for elongated embryos, 
regardless of bulk phase stability without the field. 

In the second term of Eq. ([5]), we allow fj, (included in 
Ro) to be negative for a new phase that is energetically 
unfavorable in the bulk. Eq. ([5]) predicts needle-shaped 
second phase particles to be energetically favorable pro- 
vided that, 
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Figs. [5] shows, indeed, that in a strong field the system 
lowers its energy more easily by forming elongated par- 
ticles, regardless of the sign of \x. That is the general 
mechanism by which the field can drive the transition. 

While post-nucleation growth is beyond the present 
scope, the final state can be readily described in the same 
framework. If the final state is a stand-alone metallic 
body in the same field E, then the obvious modifications 
to Eqs. (f5|) and © will be as follows: [i will have the 
meaning of the chemical potential of that metal, a will 
stand for its surface tension, and e should be set to unity 
assuming that the body is in vacuum. Therefore, the 
transformation will result in a uniformly metallic needle- 
shaped body. 

Eq. (f5]) suggests that nuclei with R — > are most 
favorable. Realistically, R must be greater than some 
minimum value determined by extraneous requirements, 
such as sufficient conductivity to support a large dipole 
energy or mechanical integrity. Based on relevant data, 
it was estimated Q that R m in — ctRo, where a ~ 0.1 is 
a phenomenological parameter. That puts R m in in the 
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range of molecular size. The free energy in the region 
R < Rmin is substantially larger than described by Eq. 
([5]) because the energy reducing effect of the electric field 
cannot be manifested by such thin particles; this can be 
approximated by a potential wall. With the latter in 
mind, the maximum of the free energy in Eq. ([5]) (with 
R = ai?o) yields the nucleation barrier 



W C yl = WO' 



E ' 



(7) 



The associated critical aspect ratio is, H c /R m i n = 
Eo/iEa 1 / 2 ) > 1. 

The barrier of Eq. Q is suppressed when E > E c = 
a 3 / 2 E . Correspondingly, nucleation of needle-shaped 
particles is vastly accelerated by electric fields under 
which spherical particle nucleation would be practically 
unaffected [cf. Eq. flHJ], as illustrated in Fig. [3l More- 
over, for the case of /i < 0, there exists a field range 
(E < Eq) where spherical nucleation is not possible and 
nucleation only occurs via needle-shaped embryos (left of 
the vertical line in Fig. [3]). 
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FIG. 3: Normalized nucleation barrier, W/Wo, as a function 
of the applied field, E, relative to the critical field E c . The 
barrier for nucleation of elongated particles (red, solid line) 
[Eq. 0] is compared to that of spheres [Eq. @] for the 
cases of: a stable new phase \p, > in Eq. Q] (blue, dashed 
line); and an unstable new phase [fj. < in Eq. {T}] (blue, 
dash-dot line). Nucleation of elongated particles is highly 
favored across the entire field range (E c < E < Ec/o?), where 
a — 0.1. Nucleation in the region E > E c /a 2 is uncertain 
due to the requirement of ultra-small nuclei. In the region to 
the left of the vertical line, the field can drive the transition 
to a phase that would be unstable without the field; that is 
the region where otherwise unobtainable materials may be 
synthesized. 

It should be noted that while significant as the phase 
change driver, the electric polarization here remains 
small with respect to the total charge distribution in 
the needle-shaped embryo. Indeed, the charge moved 
to the embryo ends is estimated as ~ EH 2 , with 



H c given below Eq. ([7j). That should be compared to 
the total charge of the embryo q ~ eH c R 2 min jc? , where 
e is the electron charge and a is the characteristic in- 
teratomic distance: qE/q ~ E$j (E at \fa) <C 1, where 
E a t ~ e/a 2 ~ 10 9 V/cm is the characteristic atomic field. 
A more accurate analysis in Ref. (p. 17) shows that 
the ratio qE I q is further reduced by a factor of 1/L, where 
L > 1 is defined in Eq. gj). 

Another comment concerns field-induced transitions in 
the vicinity of a bulk phase transition determined by 
the critical temperature (T c ), pressure, or concentration, 
such as e. g. the bulk phase transition between the insu- 
lating and conductive phases of VO2 at T c — 340 K. Us- 
ing the standard approximation fj, = fio(l — T/T c ), where 
/io is the chemical potential difference between the two 
phases at zero temperature, results in the corresponding 
renormalization [cf. Eq. ([5])]. 

W cx (1 -T/T c y 2 ,R cx (1-T/T C )-\E cx {1-T/T C ) 1 ' 2 . 

Because Rmin is determined by the microscopic struc- 
ture and remains practically independent of T, we ob- 
serve that a cx (1 — T/T c ). As a result, the barrier 
W cy i is temperature independent. That conclusion is 
in striking contrast to the prediction of CNT that the 
nucleation barrier is strongly temperature dependent, 
W cx pr 2 cx (1 - T/T c )~ 2 . Thus, we observe that field- 
induced nucleation becomes exponentially more effective 
than the classical nucleation of spherical particles in the 
proximity of bulk phase transition. It can dominate even 
under relatively weak fields E > E oaQQ 2 (l — T/T c ) 2 , 
where Eqo and cvoo are the zero-temperature values of 
Eq and a. This effect can be rather substantial. For 
example, 1 — T/T c s» 0.2 for the case of VO2 at room 
temperature yields E > 100 V/cm. We note that needle- 
shaped nuclei in polycrystalline VO2 have been observed 
[llT | , which can be attributed to nucleation in the internal 
fields induced by the grain boundaries. 

Field driven phase transitions would also be enhanced 
in laser or dc fields that are sufficiently strong to ionize 
the material. Indeed, that process would generate free 
charge carriers, thereby increasing the system polariz- 
ability and its related trend toward the transformation. 

As a provocative example, consider next the synthe- 
sis of metallic hydrogen (MH) . Predicted by Wigncr and 
Huntington [l2| in 1935, solid MH has not yet been ob- 
served under static pressures of up to 342 GPa [HQ. 
Dynamic c omp ression beyond 200 GPa has also been em- 
ployed 15|, [l6( . The only direct evidence thus far was the 
brief observation [Tj]] of a highly conductive liquid phase 
under a Shockwave pressure of 140 GPa and temperature 
around 3000 K. We will now attempt a rough estimate 
of the electric field range under which MH could be syn- 
thesized under standard pressure. 

We use fj, ~ 0.1 Ry/atom [l8| for the difference in 
chemical potential between the molecular and monatomic 
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phases, and a ~ 1 Ry/atom (as a rough order of magni- 
tude estimate for significantly different structures) , yield- 
ing R ~ 10 A and W n ~ ^i? 3 , - 10 3 eV. Assuming 
q = 0.1, and e ~ 1 [15], we obtain the critical field 
E c = c?l 2 E?> ~ 10 7 V/cm, or equivalent laser intensity 
I c ~ 10 12 W/cm 2 . Therefore, the practical window for 
field-induced synthesis of MH is 10 7 < E < 10 9 V/cm. 
That field range could be made even lower in the proxim- 
ity of the bulk phase transition (e.g. close to the critical 
pressure). Investigation could also be conducted with 
hydrogen rich alloys, such as CH4 (or other paraffins) or 
Siti4 (112)2 [HI . We note that while MH particles can be 
field-induced according to our estimates, they will remain 
metastable and will exist for only a finite time upon field 
removal due to the inequality in bulk chemical potentials 
(fx < 0) and the metastability barrier of ~ 1 eV [2(| • 

We shall end by briefly mentioning the possibility of re- 
verse, metal-to-insulator transformations when the elec- 
tric potential (rather than the electric field) is kept con- 
stant. Here we limit our arguments to an analogy with 
the well-known elementary physics problem of a capac- 
itor with plates partially immersed in water. When it 
is disconnected from the voltage source U, its charge Q 
(hence, electric induction) is conserved and the energy 
Q 2 /2C tends to decrease via the increase in capacitance 
C ; this is achieved by pulling in the more polarizable sub- 
stance (water). This is analogous to creating more polar- 
izable metal particles in the above analysis. However, the 
capacitor will push water away when it is connected to 
the voltage source in order to decrease its energy CU 2 /2. 
In terms of our consideration, this corresponds to elimi- 
nating the more polarizable metal phase. 

In conclusion, we have shown how a symmetry- 
breaking electric field can drive an insulator-metal phase 
transition in any dielectric through nucleation of needle- 
shaped particles, eventually transforming it to a uni- 
formly metallic elongated body. We determined the con- 
ditions under which such phenomena are possible and 
their corresponding transformation rates. From a prac- 
tical standpoint, the concept of field-induced phase tran- 
sitions can stand as a unique pathway for synthesis of 
materials that may be otherwise unobtainable. 
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